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Abstract: A phase of massive gravity free from pathologies can be obtained by couphng 
the metric to an additional spin-two field. We study the gravitational field produced by 
a static spherically symmetric body, by finding the exact solution that generalizes the 
Schwarzschild metric to the case of massive gravity. Besides the usual 1/r term, the main 
effects of the new spin-two field are a shift of the total mass of the body and the presence 
of a new power-like term, with sizes determined by the mass and the shape (the radius) of 
the source. These modifications, being source dependent, give rise to a dynamical violation 
of the Strong Equivalence Principle. Depending on the details of the coupling of the new 
field, the power-like term may dominate at large distances or even in the ultraviolet. The 
effect persists also when the dynamics of the extra field is decoupled. 
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1. Introduction 

The search for large-distance modified theories of gravity, motivated by the evidence for the 
cosmological acceleration, has stimulated a number of studies in the recent years. The main 
goal has been to look for a massive deformation of standard general relativity, featuring a 
large distance (infrared) modification of the Newtonian gravitational potential, and massive 
gravitons. 

The idea of considering a Lorentz-invariant theory of a massive spin-two field dates 
back to 1939 [|^: the resulting theory is plagued by a number of diseases that make it 
unphysical, besides being phenomenologically excluded. In particular, the modification of 
the Newtonian potentials is not continuous when the mass vanishes, giving a large 
correction (25%) to the light deflection from the sun that is experimentally excluded 
A possible way to circumvent the physical consequences of the discontinuity was proposed 
in the idea is that in the Fierz Pauli theory (FP) the linearized approximation breaks 
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down near the star and an improved perturbative expansion must be used, leading to a 
continuous result when m — > 0. Whether the solution associated with the improved pertur- 
bative expansion valid near the star can be extended up to infinity is an open problem Q. 
In addition, FP is problematic as an effective theory at the quantum level. Regarding FP as 
a gauge theory where the gauge symmetry is broken by a explicit mass term m, one would 
expect a cutoff A2 ~ mg~^ = (mMpi)^^'^ , however the real cutoff is A5 = (m^Mp^)^/^ ||5| 
much lower than A2. A would-be Goldstone mode is responsible for the extreme ultraviolet 
sensitivity of FP theory, that becomes totally unreliable in the absence of proper comple- 
tion. These issues cast a shadow on the the possibility of realizing a Lorentz-invariant 
theory of massive gravity |Q. 

It was recently noted that by allowing lorentz-breaking mass terms for the graviton 
the resulting theory can be physically viable being free from pathologies such as ghosts 
or low strong coupling scales, while still leading to modified gravity. Since the mass terms 
break anyway the diffeomorphisms invariance, this possibility was analized mainly in a 
model-independent way, by reintroducing the goldstone fields of the broken gauge invari- 
ance, and by studying their dynamics |5|, |8[; we refer to a recent review for the status and 
results in this direction [^. This approach has the power of being model independent, but 
the advantage turns into a difficulty when investigating the concrete behavior of solutions. 

In 1 10] we considered a class of theories that generate Lorentz-breaking mass terms 
for the graviton, by coupling the metric to an additional spin-two field. This system was 
originally introduced and analyzed by Isham, Salam and Strathdee 12] and reanalyzed 



more recently in ]|13, 14 , ^q] . While this approach may seem antieconomical, we stress that 



this is the simplest model that can explain dynamically the emergence of lorentz breaking 
and give mass to the graviton. What happens is that the two tensor fields lead in general 
to two coexisting and different backgrounds, inducing Lorentz-breaking mass terms at 
linearized order. For a general discussion on the consequences of Lorentz Breaking see |l6| . 

The linearized analysis showed that only two gravitons propagate, one massive and the 
other massless, both with two polarization states, representing two kinds of gravitational 
waves (GW). These are the only states in the theory that feel the Lorentz breaking, showing 
a frame-dependence that may be measured at future GW detectors. 

In addition, the linearized gravitational potential differs in a crucial way from the 
Newtonian one: it contains a new term that is linearly growing with distance. Of course, 
this signals the breakdown of perturbation theory at large distances, and in this regime 
the theory should be treated fully nonlinearly. This fact is not surprising, since one has 
effectively introduced nonlinear interactions, and therefore antiscreening may be present 
also at classical level like in non-abelian gauge theories. We believe that this is a general 
feature of massive gravity theories due to the presence in the full theory of nonderivative 
interaction terms. In such situations the linearized analysis is of limited reach, and we are 
forced to find exact classical solutions to be compared with the standard Schwarzschild 
metric. Though in general this a very hard task, we have managed to find a whole class of 
interaction terms for which nontrivial and rather interesting exact solutions can be found. 

After describing the setup and the flat backgrounds in section |2| ad 2.1 , 



we review 



the linearized analysis and its problems in section 2.2. We then describe the spherically 
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symmetric solutions in section |^ that we match with an interior star solution to estimate 
the modifications of the gravitational potential as a function of the source parameters. 
We also comment on the properties of these solutions and of the whole theory in the 
interesting limit when the second metric decouples, leaving just one massive gravity theory 
with modified Schwarzschild solutions, as well as in the Lorentz-invariant limit. 

2. The model 

Consider a gravity theory in which, besides our standard metric field, an additional rank-2 
tensor is introduced in in the form of a bimetric theory. The action is taken as^ 

S = j(fx {Mli^Ri + £i) + {Mli^R2 + £2) - ^{gig2f'^V{X)\ , (2.1) 

and for symmetry each rank-2 field is coupled to its own matter with the respective La- 
grangians Ci^2- In the interaction term we only consider non-derivative couplings. The 
only invariant tensor, without derivatives, that can be written out of the two metrics 
is Xy = gi°'g2av^ ^'^d then V is taken as a function of the four independent scalars 
{r„ = tr(X"), n = 1,2,3,4} made out of X. The cosmological terms can be included in 
V, e.g. Vk, = Alr^/^ with q = deiX = g^lgx. 
Then the (modified) Einstein equations read 

Mp2,iEii; + Qi(; = iri(; (2.2) 

Mp'z2^2(; + Q2i; = ^T2^ (2.3) 

where we defined the effective energy-momentum tensors induced by the interaction: 

Q^M = q^l^ ^b^y - ^{y'XX\ (2.4) 
Q2/; = g-i/4[y5/. + 4(y'j^)/;]^ (2.5) 

with {y'yt, = dVjdX^^. Indeed, the field 52 plays the role of matter in the equations of 
motion for gi, and viceversa for g\. 

The Einstein tensors satisfy the corresponding contracted Bianchi identities^ 

gT^i^E^^, = V\E^^, = gr^2aE2^u = V^i?2^. = . (2.6) 

that follows from the invariance of the respective Einstein-Hilbert terms under common 
diffeomorphisms 

5gi^, = 25i„(^Vi,)r ^52^. = 252„(^V2,)r • (2-7) 

The interaction term is also separately invariant and we can derive conservation laws for 
Qi and Q2 similar to the conservation of the energy-momentum tensor in GR: 

^iQipiv = on shell for g2 

(z.oj 

'^'iQ'i^iu — ^ shell for gi . 



^We use the mostly plus convention for the metric. Indices of type 1(2) are raised/lowered with gi{g2)- 
denotes the covariant derivative associated to the Levi-Civita connection of gi/2- 
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These identities are quite powerful; for instance they ahow to solve completely the sim- 
plest of these models, when V is a function of q only. This peculiar case is discussed in 
appendix 

2.1 Asymptotic solutions 

At infinity, far from all the sources, we expect that gi and g2 are maximally symmetric, 
setting up the benchmark for the asymptotic behavior of all solutions of the EOM. Denoting 
with — /Ca/4 the constant scalar curvature of gafiu, i-e. 

M^laEaf^u = ICa Qa^^u , = 1, 2) (2.9) 



the equations (|2.2| )-(2.3) read 

2V + {q-^'^lCi + gi/^/Ca) = (2.10) 
8 (V'XY + 51^ (g^/^/C2 - q-^'^Ki) = . (2.11) 



and these equations can be solved for specific ansatze. In order to study the properties of 
this model for asymptotically flat spaces, we analyze first the biflat solutions, /Ci = /C2 = 0. 
Eqs. (|1D-(|T1|) yield: 



y'^ = 0, y = 0. (2.12) 

Assuming that rotational symmetry is preserved and that the two metrics have the same 
signature, the biflat solution can written in the following form 

91 = VfMu = diag(-l, 1,1,1) 

(2.13) 

92,iiy =^'^ diag(-c^, 1,1,1), 

where c parametrizes the speed of light in sector 2 and lv is the relative conformal factor. 
Eqs. ( p. 12 ) correspond to three independent equations y = 0, V'^ = and V'l = 0, where 



and i = 1,2,3 stand for temporal and spatial indices. Therefore, two of these equations 
determine the values of two parameters c and lo, while the third represents a fine tuning 
condition for the function V, necessary to ensure flatness. The same sort of fine tuning is 
necessary in the context of normal GR to set the cosmological term to zero. Therefore, for a 
generic function V{X) we expect to have a Lorentz Breaking (LB) solution with c 7^ 1, and 
hence a preferred reference frame ( |2.13| ) in which both metrics are diagonal. Asymptotic 
flat solutions should approach ( |2.13| ) in a suitable coordinate system. In addition, there 
exists also a Lorentz Invariant (LI) solution, with c = 1: in this case two equations coincide 
V'q = V'l, and can be used to determine the value of lo. 

Summarizing, asymptotically the solutions fall in two branches: LI with c = 1, and LB 
with c^\? The LB branch is of particular interest since it naturally allows for consistent 
massive deformations of gravity p 



^In the special case when V — l^(detX), Bianchi identities force detX to be constant, and there are 
additional gauge symmetries that allow to set c = 1. As a result the branches are equivalent (appendix 
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2.2 Review (and critics) of the linearized analysis 



In [10 1 we performed a linearized analysis around the biflat background, that we report 
here for the LB branch. 

In addition to the kinetic terms, the linearized action contains a Lorentz-breaking mass 
term for the fluctuations hai_iv = Qa^iu — da^iu (o = 1)2). Since on the biflat background 
V = V"l^ = 0, one can expand the potential at second order in the fluctuations Y = 
^g2^^2 ~ di^hi X, and define the mass lagrangian: 

jCm = -2(9192^^ tr [YV"{X)Y] 

= \ (^00 -^0 /ioo + 2/ioi Mx hoi - h\j M2 hj + 4 -^^3 ha - 2hi M4 /iqo) • (2.14) 

Here h^i, = {hi^^, h2f_iu} is the column vector of fluctuations and A^o,i, 2,3,4 are 2x2 mass 
matrices. It is then crucial to realize that, due to linearized gauge invariance (that we 
remark is never broken) these matrices are of rank-one; one can write 

Mo = XoC-^VC-^ 

M2,3 = A2,3 r, ^ ^ "1^) ' ^ ^ ( ^ C ) ^^-^^^ 

where Ao,2,3,4 depend on the potential. In addition, due to the LB, M.i vanishes regardless 
of the potential V . This fact leads to a well defined phase of linearized massive gravity.'^ 

In this phase, the only propagating states are the two spin-2 tensor components of 
the fluctuations (two polarizations each) corresponding to two gravitons, of which one is 
massless and the other has mass A2. Their dispersion relation is non-linear due to their 



mixing and their different propagating speeds [10|. 

The other components, scalars and vectors, do not propagate, and therefore discon- 
tinuity and strong coupling problems are absent in this phase. They however mediate 
instantaneous interactions, so the Newtonian potentials that one finds at linearized level 
are then drastically modified; for example in sector 1, the potential from a point-like source 
Ml is: 

^^ = -9Ml^GMi^?r, (2.16) 
r 

where 

2 _ A2 3A4 - Ao(3A3 - A2) (o^7\ 
^ -2Ml Al-Ao(A3-A2) ^ 

and for later reference note that may be negative. 

The linearly growing term in ( p. 161) signals the breakdown of perturbation theory at 
distances larger than r//j = {GMiy?)^^ ||8|, |l^], and one usually considers the solution to 
be valid as long as the potential stays in the weak field regime. However, one should note 



*The vanishing of the second eigenvalue of A^i can be understood by noting that hioi~h20i is a goldstone 



direction, corresponding to the broken boosts in the LB background. In sec. 3.6 we comment on the fate of 
this condition on nontrivial backgrounds. 
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that the hnear term in ( p.l6 ) is induced by an other scalar field having an instantaneous 



interaction and acting as a source for <I> (see e.g. |22, ^). It is then easy to realize that non- 
linear corrections to this field can drastically modify the IR behavior, even in the weak-field 
regime. We can clarify this point by showing, as an example, two systems of differential 
equations that differ by nonlinear terms and have drastically different IR behavior 

f A$ + /iV = M6^ (x) f A$ + /i^fj = M6^ (x) 

[Aa = M6^{x) [Aa + Xa"^ = M6^{x), ^' ' 

Here $ is a scalar field (mimicking the gravitational potential) a is an additional scalar 
field coupled to it by a mass term (and A is the laplacian). While in the first system 
a ~ M/r and this induces a linear term in <I> like in (2.16), in the second system a drops to 



zero faster than M/r so that the bad behavior of ^> is cured. What happens is that the IR 
behavior is dominated by a non-linear term, because effectively A ^ at large distances. 
We incidentally point out that standard GR is safe in this respect, because nonlinear terms 
coming from the Einstein tensor are always accompanied by two derivatives and thus are 
equally suppressed at large distances. 

We are thus led to the conclusion that in massive gravity the situation is similar to 
non-abelian gauge theories, where the large distance behavior is generically non-trivial and 
inaccessible to the linearized approximation. We recall that in Yang-Mills theories, non- 
abelian configurations of charges can lead to non coulomb-like classical solutions, screening 
or even anti-screening the charge, leading also to infinite energy configurations p!7[| . 

In this situation what one may try is to really look at higher orders and maybe retain 
the first terms that are relevant at large distance. This approach would require the painful 
procedure of defining the gauge invariant fields at higher orders, and would also lead to 
nonlinear terms mixing scalars, vectors and tensors. Instead of following this approach, we 
find more instructive to study the exact spherically symmetric solutions. 



3. Exact spherically symmetric static solutions 

The Schwarzschild solution describes the spherically symmetric gravitational field produced 
by a spherically symmetric source. It is crucial to understand what kind of modification is 
introduced in this theory by the presence of a new spin 2 field. Spherical symmetry allows 
us to choose a coordinate patch {t,r,9,{p) where gi and 52 have the form 

dsl = -J df + K dr'^ + r'^ dn'^ (3.1) 
dsl = -C dt^ + Adr'^ + 2D dtdr + B dn^ . (3.2) 

and all the functions J, K,C, A, D, B entering gi and 52 are function of r only. Notice that 
the off-diagonal piece D cannot be gauged away. 

3.1 Black hole solutions 

In the absence of matter, a number of interesting properties follow from the form of the 
Einstein tensors Ei^, derived from ( |3.1| )-( |3l2|) and do not depend on the chosen V. 
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Following |T^, 15] the spherically symmetric solutions can be divided in two classes: type 
I with D 7^ and type II with D = 0. We shall focus here mainly on type I solutions. 

Since is diagonal by the choice of the first metric, then also {V'X)p must be 
diagonal because of the EOM (p.2|). The only possible source of a off-diagonal term in the 
RHS of ( ^ ) would be (V'X)'^, so as a result also E2u must be diagonal, i.e. £'2* = 0- For 
type I solutions, this condition amounts to a single equation: 

AC + D^ = d2^-^, (3.3) 

where d2 is a constant. Incidentally using this relation it turns out that E2I = £2^, then 
using Q also (y'X)l = (^'X);, and by Q we have also that Eil = Ei^. This relation 



determines K in terms of J 

= f , (3.4) 

with di an other constant. 

The metric 2 can be brought in a diagonal form by a coordinate change dt = dt' + 
drD/C. Thanks to ( p. 31 ), in the new coordinates we have 

dsl = -Cdt'^ + ^^dr^ + Bdn^ (3.5) 
4B C 

(and of course the metric 1 in no longer diagonal). Then by a suitable change of r the 
metric 2 can also be put in a Schwarzschild-like form; setting r' = \J B{r), we find 

dsl = -C dt'^ + ^ dr'^ + r'^ dn\ C{r) = C{r'), (3.6) 

which shows that C is the physically relevant potential in sector 2. 

To proceed further a choice of V is needed. In the existing literature essentially all the 



results are based on a potential Vis introduced in and 1 11 ] in the context of hadronic 
physics.^ The motivation for this choice is probably due to the fact that Vis is the simplest 
potential producing a FP mass term in the (Lorentz-invariant) linearized limit: 

Vis = {r-2 - t!i + 6t_i - 12) 

= (/^" - g'nig""' - 9'ni9i,p9iua - gi,ugipa) (3.7) 

— tr(/i^) — tr(/i„)^ for gi = g2 = V ■ 



where = /i2^i/ — hi^j^^. For Vis it was shown in |12] that type-I solutions are always 

Schwarzschild-(A)dS, and it was recently realized that these solutions are present for any 
potential [|2^. It turns out that Vis can be deformed and there exists a whole family of 
potentials for which the exact spherically symmetric solutions can be found. 
Let us consider the family of potentials 

V = ao + aiVi + 02^2 + 03 V3 + 04^1 + hV-i + b2V-2 + 63VI3 + 64^-4 

+ q-^/^Ai + q^/^A2 , (3.8) 



^In the years preceding QCD, the proposal of Isham, Salam and Strathdee was that of a second metric 
mediating a strongly coupled interaction, responsible for confinement of quarks inside tiny black holes. 
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where we introduced the following combinations involving the generalized determinants 
(again r„ = tr(X") and e is the 4-index antisymmetric symbol) 

^0 = o/if I (g€g2g2ff2g2) = 1 = ~ + 8x1X3 + 3 t| - 6X4) 

24 1 52 1 24g 

1 1 o 

^1 = ^1 — r(ee925'2S'25i) = ("T"-!) = ^-l"^! - 3T2ri + 2T3) 

V2 = -\ — I (ee925'25'i5i) = ("^^i - ^-2) = ^"^'^1 - ^2) (3.9) 

V3 = ■r^(eeg25'i5'i5'i) = (^^i - 3r_2'r_i + 2t_3) = Qq'^ n 

V4 = — Aeegigigigi) = (rl^ - 6t„2T^i + 8r_iT_3 + 3r^2 - 6r_4) = 24 g"^ 
152 1 

and where y_„ = 14(X ^ X^"^)- The cosmological constants Ai and A2 have been added 
to simplify the asymptotic flatness conditions. The Isham-Storey potential is recovered by 
setting ao = —12, oi = 6, 02 = 1 and 61 = 62 = &3 = &4 = ^3 = 04 = 0. 

Remarkably, the general combination V of (|3.8|) leads to solvable equations for type I 
spherically symmetric solutions, and these can be found in a closed form (these equations 
are the main result of the paper) :^ 



J 

C = c^u? 



1 O*^™-! , V 2 

1 — 2 h /Cir 

r 



Gm2 2 



+ 2G5rT, KJ=1, (3.10) 

2G 



5rT, Z)2 + AC = cV (3.11) 



B = u?r\ A = ^^L_^^^l!Z^^ (3.12) 

with {J,C} = {J,C}/a;4(c2 + 1), 5 = 5/A2 [(c^ - 1)(7 + 1)(7 - 2) + 1)]. 
The solution depends on the integration constants mi, 1112 and S and we have introduced 
G = I/WttM^i^ and k = M^i^/M^ii- The values of c^, 7 and of the graviton mass A2/Mp;^ 
of the linearized analysis ( |2.15| ) are given in terms of the coupling constants: 

^2 _ ai + 452 + 603 4[(a2 + 3a3)-c^(fe2 + 3fe3)] _ 2(7 - 2) ^ 

61 + 452 + 663' ^ 02(61+462+663) ' ^ 7 asj, 

where a„, = uj^'^'^an, bn = io'^"bn and a„ = (a„ — c^bn){c^ — l)/c^. Notice that one may 
trade ai, 61 and e.g. 02 for, respectively, c^, 7 and the graviton mass A2, showing that these 
may take any value for this class of potentials. When 7 < 2, the /Cj are proportional to 
the constant asymptotic curvatures of gf, the explicit expressions are given in appendix^. 
Finally, cj^ is also in general a free parameter that determines for example A2 to have 
IC2 = 0, after having fine tuned Ai to set /Ci = 0. 



®This solvability is linked to the fact that the combinations Vn are actually the coefficients of the secular 
equation of X and are (multi)linear combinations of its eigenvalues Aii Vn — X^ij>i2 - >i ^ii^i2 ' ' ' ^in- 
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The expression ( ^.lOD resembles the Schwarzschild-dS(AdS) solution but with a crucial 
difference: a r'*' term of magnitude S is present and it may alter significantly the behavior 
of the gravitational field, depending on whether 7<— 1, — l<7<0or7>0. 

Before discussing these solutions, let us comment on the Isham-Storey potential Vis 
used traditionally. Since in this case all 6„ vanish, it leads to a singular situation where 
^ oo unless an additional fine tuning u;^ = 2/3 is performed. Even choosing this case, 
the linearized analysis is ill defined due to an enhanced gauge invariance (see [10| for the 



case of A,; = 0), and moreover from ( 3.13| ) one has 7 — > oo. This is the reason why only 



standard Schwarzschild-(A)dS solutions were found. 

Now, in order to shed light on the physical meaning of the various constants in the 
solution let us also compute the total gravitational energy, as measured with respect to 
backgrounds 1 or 2. In the stationary case this is the Komar energy that can be calculated 
as a surface integral on a sphere of large radius router 00 (see appendix We find, for 
the two fields: 



£1 = mi + 5'7r2'+ter 
c 

UJ 



£2 = m2-—S-f r2^i^, . 



(3.14) 



From these expressions we see that only IR modifications with 7 < — 1 will lead to finite 
total energy. 

Case 7 < — 1: At very large distance the solution reduce to a maximally symmetric 
solution parametrized by ICi. In particular one can set /Ci = /C2 = with a single fine 
tuning, determining the asymptotic conformal factor as discussed above, so that the 
solution describes asymptotically fiat metrics. Clearly because 7 < —1, at large distances 
gravity is Newtonian, while at short /intermediate distances, depending on S, the presence 
of the additional spin 2 field has changed the nature of the gravitational force. 

Since the large distance behavior is Newtonian, the total energy is finite; taking 
fouter — > 00, we find £1 = 772-1, £2 = "i2. In black hole solutions like these, mi and 
m2 are just parameters, that can be related to the mass of a material object only when the 
solution is considered as the outer part of, for instance, a star. In the case of standard GR, 
for a star of radius R and mass density p, the total gravitational energy E is the total mass 
M = 4ttR^p/3. Here, the interaction with g2 is turned on and we expect a contribution to 
this energy given by the interaction term Qi. Its size should be controlled by V and by the 
matter itself, because this interaction energy also is turned on by the source. Moreover, 
by dimensional analysis the coefficient S of the r"^ term should also be a function of the 
size of the object, and not only on its mass. This can be understood intuitively as the 
failure of the Gauss theorem due to the presence of Q in the EDMs and of the term in 
the solution. Accordingly, the separate contribution of Qi to the energy is not expressible 
as a fiux on a 2-surface at infinity, as it happens for the total Komar energy. The explicit 



computation of this interaction energy for a star will be performed in section 3.3 



Case 7 > — 1: For simplicity, also in this case we set /Ci = /C2 = as discussed above, 
but note that because the new term induces a curvature R ~ r'^~'^, only when 7 < 2 we 
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have that gi and 52 are asymptoticaUy flat and can be interpreted as an asymptotical 
conformal factor. For these choices of 7, we have a solution such that Qi does not vanish 
rapidly as r ^ 00, and compensates a slow fall-off (or rise!) of the gravitational field. 
However on dimensional grounds any fall-off slower than 1/r makes the Komar total energy 
infinite, and indeed when router 00 both £1 and £2 diverge making this configuration 
physically unfeasible. 

If spherically symmetric solutions of infinite energy are surely not physical, this may 
only suggest that solutions will not be spherically symmetric, as it happens in non-abelian 
gauge theories. For example one may speculate that finite energy configurations will ar- 
range in flux tubes of gravitation at large distance, between sources of type 1 and 2, as 
suggested by the different signs of 5" in J and C. In a similar 'confinement-like' scenario, the 
term r"' may be screened dynamically by the self-arrangement of configurations of matter 
1 and 2, so that effectively 5 ^ at large distances, as suggested by the full star solutions 
that we will describe later. 

To summarize, we found that exact black-hole solutions are modified in the IR or in 
the UV depending on the choice of the potential, and that this behavior is not captured 
by the linearized approximation. There are even cases where the behavior of the potential 
is not modified at all with respect to GR (e.g. 7 = 2 or 7 = —1) while the linearized 
approximation still shows a linear term. 

It is also interesting to observe that in the limit in which the second metric decouples, 
Mpi2 — > 00 (i.e. K — > 00), and assuming that 7712 and S remains finite, from the solu- 
tion ( |3.10D we find that the term S r"^ remains in gi while (72 becomes exactly flat. We will 
discuss the decoupling limit as well as the limit ^ 1 in section 3.4 and 3.5. 



3.2 Comparison with the linearized solution 

It is interesting to comment on the perturbative origin of the exact solutions. This can be 
addressed by looking at the asymptotical weak- field limit of the solution (for 7 < —1): 

J, K, C, Ar^ const + 0{l/r) D ^ 0{^/lp) . (3.15) 

The crucial observation is that D vanishes more slowly (and non-analytically) than the 
other components of the perturbations. As a result, this solution is not captured by the 
standard linearization, where all the perturbations have the same large distance fall-off. 

Technically, the origin of this behavior can be traced back to the equation (Q1.2)* = 0, 
that is algebraic: 



D{r) 



A{r)C{r) + L>(r)2 ai + Aa2r'^B{r)-^ + Qa^r'^B{r''-'^ 



+ 



J{r)K{r) bi + Ab2r~^B{r) + Qh^r-'^B{r) 



(3.16) 



This equation can be solved either with D = (type-II solutions) or with D ^ Q (type-I). 
In this last case, for the exterior solution, since it turns out that B = uj'^r'^, asymptotically 
the equation turns into the definition of the speed of light of 52 (as in ( 3.13D ), while the 
deviations give the mentioned behavior of D ~ 1/ y/r. 
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Prom equation ( 3.16| ) we can also understand that standard linearization (around the 
LB background) can not distinguish between type-I and type-II at leading order: consid- 
ering the standard perturbative expansion (with parameter e) where D ~ e, this equation 
starts from order e^. Also, at first order D can be gauged away: it does not appear neither 
in the linearized Einstein tensors (due to separate gauge invariances) nor in the mass terms 
(due to A^i = 0). At higher orders however one must choose = 0, otherwise there is 
a constraint on the fields A, C, J, K, B that are already determined at previous-orders. 
We reach the conclusion that the standard perturbation theory around the LB background 
may only approximate the solutions in the type-II branch (if any exist: we recall that 
nontrivial type-II solutions are not known). 

On the other hand, it is interesting that the r"' term can be recovered in a semi- 
linearized approach, where one solves exactly equation ( |3.16| ) and treats the remaining ones 
perturbatively. This will be done for the interior star solution and the result containing 
the r"^ terms can be found in appendix]^, e.g. equation B.2. Alternatively, if one insists 
in solving perturbatively all the equations, the correct result can also be recovered by 
assuming D ~ -^e while all the other fluctuations are still of order e, and retaining the first 
nonvanishing order.'' 

Exactly as in the comparison between the Newtonian and Schwarzschild solutions, the 
final difference between this semi-linearized and the exact solution is just that K = 1 + 2^ 
instead of K = = 1/(1 — 2$) (and similarly for §2). 



3.3 Interior solution 

In order to determine the integration constants mi,m2 and S one can imagine that ( 3.1C| ) 
is the exterior portion of the solution describing a spherically symmetric star. We aim at 
finding the interior solution and then determine mi, 7712, S by matching with the exterior 
one. 

It is instructive to consider first in full generality a spherical star made of fluids of 
type 1 and 2, extending from the origin to radii Ri, R2, stationary with respect to the 
respective metrics: 



/- 



V 



■pi 



\ 



Pi 



Pi 



(-P2 %{P2+ P2) 



nr. 



pij 



P2 



P2 



\ 



\ 



P2 I 



(3.17) 



Like in the vacuum, since g\, T\ and E\ are diagonal, so should be Q1/2, i-e. Qi* = 0. 
This equation, being the same as in the vacuum case, is exactly solvable for the class of 
potentials (|3.8| ). The remaining equations are more involved in the presence of matter, but 
in linearized approximation the solution can be found analytically. 



similar approach was envisaged in to find an asymptotically flat modified Schwarzschild solution 
for LI massive gravity, valid in the limit. Also in that case a field that is not determined at linearized 

level for m = 0, is found to vanish non-analitically as sjljr. However that solution is not valid beyond 
some distance scale, and there is probably no global extension j^. In the present work, it is remarkable 
that the semi-linearized solution is also extendable to the exact one. 



- 11 - 



According to the discussion of section 3.2 this partial hnearization corresponds to 



choosing the type-I class {D ^ 0) also for the interior solution. 

For simplicity, we consider a star made of an incompressible fluid of constant density 
and small pressure, p <C p. The interior solution is then matched by requiring continuity 
of C, J, B, K and of the derivatives C", B' . This procedure, in the physical case when 
7 < — 1, determines exactly the exterior constants mi, m2, S.^ 

While the detailed solution is given in appendix]^, we present here the instructive case 
Ri = R2 = R, Mi,M2 7^ 0, and then discuss the phenomenologically interesting case of 
only matter 1, M2 = 0. 

Case with both kinds of matter: For Ri = R2 = R, and setting Mi 2 = 47r/9i 2^?! 2/3, 
the matching condition gives: 



mi = Ml + AM , AM = a /i^i?^ Mi M2 

m2 = M2 - AM/cKW^ (3.18) 
S= AMi?-(T+i) 15/(27 -1)(7- 4), 
where a = Sc^^'^tu'^ /5{'j + 1)(7 — 2) and /x^ is the same constant that appears in the 



linearly growing potential (2.16) of the linearized analysis. In the exterior solution ( 3.1C| ) 



S and AM modify the form of standard Schwarzschild solution. The first modification 
is in the Newtonian terms, and amounts to a mass shift with respect to the standard 
values mi^2 = ^1,2- The second is the new term r'^ . Both are proportional to the same 
combination AM.^ 

The mass shift AM can be understood as the contribution of the interaction terms 
Q to the total energy, i.e. to the total mass as measured by the Newton law at large 
distance. To clarify this, it is useful to recall that in standard GR the Komar energy, 
written as a spatial volume integral (Svr)^^ J R^dv (see appendix ^ can be rewritten 
as a volume integral of the matter energy-momentum tensor by means of the Einstein 
equations: E = £t = (87r)-i/(2r^ - T5'p|i^'dvu. This resuh is modified m massive 
gravity because the Einstein equations contain the additional energy-momentum tensor of 
interactions Q, and the additional contribution can be evaluated with its volume integral 
£q = (Stt)"^ fi'^Q'^ ~ Q'^p^'^dvj,. We remark that while the sum of these two integrals, 
being the total energy, can be expressed as a surface integral at infinity (see appendix 
they separately can not, and they can only be evaluated using the smooth interior and 
exterior solution. The result of the volume integral is (again finite only for 7 < —1): 

= = ^1,2 , = -k£q, = AM . (3.19) 

This confirms that the mass shift AM is a screening effect, due to the energy of the 
interacting fields in Q, and corresponding to the nonzero Ricci curvature even outside the 
source (see e.g. [^]). 

*In the unphysical case 7 > — 1, although the potentials J and C are regular, the field B develops a 
singularity r'^'^^ in the origin, calling probably for a fully nonlinear interior solution. 
^But see appendix ^ for the full case _Ri ^ R2. 
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As a side remark, looking at the matching ( |3.18| ), we observe that we did not hnearize 
in V , but neglecting terms higher order in the matter density one has effectively neglected 
higher orders in y ~ Indeed, the result depends on the two dimensionless parameters 
E?^'^ and GM/R^ but at first order in GM/R only the first order in R?^"^ appears, i.e. 
GMRfj?, and the final result is smooth when the interaction vanishes, 1/ ~ /i^ ^ 0. This 
is opposed to the singular massless limit of Lorentz-Invariant (Fierz-Pauli) massive gravity. 

Case of normal matter: Turning off M2 = 0, we can focus on sector 1 and discuss 
more phenomenologically how normal gravity is modified by the presence of the additional 
spin two field. 

From the matching condition we have, with Mi = M: 

mi =M{1 + af/R^) , m2 = -a R^ M /c ku? 

(3.20) 

S=fi^MR^^^ 15a/(27-l)(7-4) 



and from ( p. 10 ) we find for the modified potential (ignoring the numerical factors): 



kl + ^^^R^) + ^^R{^y 



(3.21) 



The mass shift is now equivalent to a rescaling of the Newton constant G{1 + afj? R'^), that 
depends on the source radius!! 

We observe that for the sun^'' we have fi^R'^ ~ lO"^*', assuming all coupling constants 
to be of the same order so that fj, ~ rUg < (lO"^'^ eV) ~ (lOOAU)^^ (this limit corresponding 
to the rough experimental bound on the graviton mass from pulsar GW emission [|8|, |2^, |l^]). 
We thus see that for the sun the size dependence is negligible and unobservable, and even 
more so for the planets. The effect becomes important for objects of size fi~^. For instance, 
for large objects with R > 10^ Rsun (red giants, large gas clouds, galaxies. . . ) the effect 
may be of order one, and induces a macroscopic modification of the Newton constant. For 
low density objects that we consider here, this modification does not depend on the mass 
but just on the object size; therefore, given the mass, a large sphere of gas has a larger 
effective newton constant. In the limit fx^R^ > 1, the surface potential would even scale as 
R'^, instead of R^ as in standard gravity. Moreover, remembering that //^ may be negative, 
the negative interaction energy could cause large fluids to antigravitate, hinting toward the 
acceleration of the cosmological solutions. 

Then, the new term in the potential is of the form 

6<l> ^ GM fi^R(^^y , (3.22) 

replacing the linear term G M fx^r of the linearized analysis. The Newtonian and the new 
term will be competing at a critical distance Tc that also depends on 7: 



R 



fi^R^ 



1 + /i2i?2 



7+1 

(3.23) 



^°i?0 ~ 5 ■ 10^ Km = 5 ■ 10^^ eV"^ and Mq = 10^^ eV. 
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Of course since 7 < — 1 the relevant modification is ultraviolet, and is evident for r < rc 
(while for 7 > —1 it would be infrared, for r > rc). To estimate we observe that since 
the exponent —1/(7 + 1) is positive, one always has < R for /^^ > 0, and the critical 
distance is inside the star. This does not mean that there will be no observable effects, 
since even subleading modifications to the newton potential may be measured (for example 
modifications of the gravitational potential of relative magnitude 10~^/~^ are at the level 
the current solar-system tests). On the other hand for negative iJ? , and in particular for 
fj^B? < —1/2, one has rc > R so that in a UV region near the source the gravitational 
potential has stronger fall-off. For ^^R? ~ —1 we even find that rc becomes infinite, so 
that the region of UV modification expands to larger and larger distances! 
We can summarize the results in the physical phase 7 < — 1: 

• For sources of dimension R < /i^^, the effects are: a mass shift equivalent to a 
small Newton-constant renormalization (1 -|- fj?R'^), and a subleading correction to 
the Newtonian potential, 6^ < (r/R)'^. 

• For large sources, of dimension R > fi^^, the mass shift is more pronounced, and for 
negative fj? even the new r"^ term can become dominant in a region near the source. 

As we see, even discarding the nonphysical and possibly confining branch 7 > — 1, the 
phenomenology of these modified static solutions appears to be quite rich, and deserves a 
separate analysis to confront their features with real physical systems, e.g. modified galactic 
gravitational field, gravitation of large sources, post-Newtonian analysis. 

3.4 Decoupling the second metric 

The idea of introducing a second metric and considering its decoupling limit, to have a 
second background at hand while disposing of its fluctuations, is not new and was indeed 
considered to tackle the problem of the nonlinear continuation of the FP massive gravity Q . 
However, due to the singular Isham-Storey potential, or due to the ill-defined nature of the 
Lorentz-Invariant theory, this did not lead to significant advance. In this work we found 
some nonperturbative solutions of the full system, so we are in a position to control the 
decoupling limit Mpi2 00 {k ^ 00) in which the second gravity is effectively switched off. 

First, as far as the propagating states are concerned, we recall from the linearized 
analysis [l^] that in the flat background out of two gravitons only the first graviton survives 
the decoupling limit: it is massive (with two polarization states and mass GX2) and has a 
normal dispersion relation. 

For the nontrivial solutions, as anticipated, one may take this limit in the exterior 
solutions, once one checks that mi, 7712 and S stay finite. This is indeed shown by the 



interior solution ( 3.18 ), therefore we directly find the result: 



= -^{1 + R'fi'a) + GMl^i'R ^-j [15a/(27 - 1)(7 - 4)] , (3.24) 
$2 = 0. (3.25) 

Both the mass shift and the new term remain, but the second gravity disappeared: here 
the other metric is flat! 
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A look at the exact solution (|3.10 ) in the decoupling limit shows that the limiting 
metric 2 is still nondiagonal {D ^ 0). This means that (72 is only gauge-equivalent to rj2 = 
w^diagj— c^, 1, 1, 1}, and that to make contact with this traditional minkowski diagonal 



vacuum one has to choose the gauge ( |3.5D , where gi is not diagonal. Explicitly we find: 

dsl = -Jdt'^ + 2Ddtdr + Kdr^ + r^dn^ , dsl = uj'^{-c^dt^ + dr"^ + r'^dO^) , (3.26) 

with J = J, K = J-^(l - D'^), D = -{cujy^J\/uj'^ - A. Notice that D is defined by the 
deviation of A from w^, and that still JK + Z)^ = 1. 

We therefore note that, to recover the present solutions in effective massive gravity 
theories, where only gi is dynamical and the Lorentz breaking is an external diagonal 
metric, one should look for nondiagonal configurations. 

We also remark that while taking the decoupling limit has left us with a flat auxiliary 
metric, still there is curvature for metric 1 in the vacuum outside the sources, due to Qi and 
Q2 being nonzero there, because of the r'^ term. Therefore the order of the limit matters, 
and one would not get the correct result if one were to assume a flat second metric before 
taking the decoupling limit. In other words, setting g2 to be flat in advance: g2 = 92, we 
have in vacuum 



El^^ = = Vb^ + 4(y'X)(; so that (3.27) 
Mli^E^^ = 2 V5^ — >V = const by Bianchi — > (Anti)deSitter . 

Instead, in the limit Mpi2 00 we have still g2f^u — > g2iiu, but different solutions for gi: 

E^^ but V6^ + 4(y'X)(; / , and then (3.28) 
^lliKv = i^^u - ^K'') + const — > non-trivial solutions. 

Summarizing, the decoupling limit shows that the theory remains well behaved, con- 
sisting of a modified gravity with massive gravitons, while the auxiliary metric is flat and 
decoupled. 

3.5 Lorentz-Invariant limit 

The Lorentz Invariant limit — > 1 is also interesting to address the Vainshtein's claim 
that nonlinear corrections actually cure the discontinuity problem in Pauli-Fierz theory Q . 
Indeed, we find that the limit ^ 1 is well behaved, and the solutions retain their validity. 
In this limiting phase therefore, gravity is modified, but lorentz breaking disappears. 

The linearized mass term is accordingly of the form a /i^^ + bh'^, however the limiting 
theory reached in this way is not the Fierz-Pauli one, where a + 6 = (and for this reason 
FP is free from coupled ghosts). Here, we get to a theory where a = 0; in fact, since 
we approach the LI phase from the Ai = branch, and because in the LI limit one has 
Xi = X2 = a, we see that also the graviton mass vanishes in this limit, as can be checked 
with the expression ( 3.13| ). It is nevertheless worth to point out that also a = is a ghost- 



free theory like a + b = 0. This case is not usually considered because at the linearized 
level there is no massive graviton as a consequence of an additional gauge symmetry (three 
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transverse diffeomorphisms) , see [23] and also the PFO phase in |1C]. Accordingly, no strong 
coupling problems are expected and no Vainshtein issues. This matches nicely with our 
model having good properties along all the LB branch, that survive also in the LI limit. 



3.6 Local Lorentz-breaking in nontrivial background 

While the asymptotic biflat metrics are Lorentz breaking, one may ask about the situation 
at finite distance. This will have definite interest when addressing the nonpropagation 
of ghosts in the described nontrivial background. In fact, we recall (section \i.2\) that on 
flat background this is a consequence of A^i = , and this follows from gauge invariance 
together with the fact that locally boosts are spontaneously broken. Now, even in nontrivial 
background a spontaneous breaking of Lorentz will lead to flat directions of the potential. 
Whether this fact will be enough to lead to absence of ghosts and to stable configurations 
is under scrutiny, and goes beyond the scope of the present work. 

To describe the local breaking of Lorentz at any given point in the nontrivial back- 
ground, one chooses a local Lorentz frame {gi = rf) and simultaneously diagonalizes g2- 
The Lorentz breaking is given by the entries of 52 j that are actually the eigenvalues of X. 
These are easily calculated (in polar coordinates t, r, 0, (p): 

X = u;^c'r\tlA}, (3.29) 

/± = 1 ± _ 



with e = i CV- + /+ + a/(c2 - l)(c2/2 - fl) 



2 

Quite remarkably, that they do not depend on the masses mi^2 of the newtonian terms, and 
this is due to the nondiagonal structure given by D. One can easily check that for r — > 00 
we have X = u}'^{(?, 1, 1, 1}, reproducing the asymptotical lorentz breaking (7 < —1). 

Then we see that in the case 5 = the eigenvalues are constant, so that at any distance 
the Lorentz breaking is the same: X = uj^{c^, 1, 1, 1}. We have two pure Schwarzschild 
solutions in a configuration that at any point breaks local boosts but preserves rotations 
(these are the solutions found in p^). 

On the other hand, since in general 5* 7^ 0, a star solution will break not only boosts 
but also local rotations, because the rr term is different from the 09 and (j)(f) ones. For 
example at large but finite distance, where Sr^~'^ is small, we have: 

X ~ uj^{c\l - Sr^-^), 1 + ^ ^ (3 30) 

To compare the situation with standard GR, we recall that in GR Lorentz-invariance at 
any given point is always valid, in the Lorentz frame, and is broken in a finite neighbourhood 
only by the curvature (tidal) effects. Here on the contrary in the gravitational sector a 
Lorentz breaking is felt also locally, and for S ^ also rotations are broken. The physical 
effect is that gravitons will propagate differently in direction of the source. 

We strongly believe that this breaking of (local) boosts and rotations at finite distance 
from a source is a general feature of nontrivial solutions in massive gravity, due to the 
presence of additional fields that can not be 'gauged away'. 
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4. Conclusions 



In this paper we approached the problem of finding a consistent massive deformation of 
gravity by introducing an additional spin 2 field g2 coupling non-derivatively to the standard 
metric field. This allows us to explore both the Lorentz invariant (LI) and Lorentz breaking 
(LB) phases working with consistent and dynamically determined backgrounds. Preserving 
diffeomorphisms and breaking Lorentz is also important; at the linearized level it forces 
A^i to vanish and no dangerous scalar mode is propagating. Still at the linearized level, 
it was shown that in the case A^i = the vDVZ discontinuity is absent, but a new 
linearly growing term is present in the static gravitational potential |1C, that seems to 
invalidate perturbation theory beyond some distance scale. To address this and the vDVZ 
discontinuity problem, we thus studied the exact spherically symmetric configurations. The 
exact solution that we found, valid for a large class of interaction potentials, shows that 
the linear term is replaced in the full solution by a power-like term r'^ , with 7 depending 
on the nonlinear couplings in the interaction potential. 

Phenomenologically, when 7 < — 1 the total energy of the solution is finite and the 
space is asymptotically flat; Lorentz is broken in the gravitational sector by the asymptotic 
value of (72) but normal matter only feels the modification of the gravitational potential. 
Using the full solution one can check that the absence of the vDVZ discontinuity is an exact 
result. The effect of the interaction manifests in the term whose size S was determined 
for a star by matching the exterior solution with an interior one. In addition to this, by the 
presence of the additional spin 2 field, the total mass of the star appearing in the Newton 
term gets a finite renormalization that depends on the object size, and may screen or even 
antiscreen the star mass. We believe that this is a general feature of massive gravity. 

When gi describes a black hole the solution depends not only on the collapsed mass 
but also on an other constant, probably remnant of the original shape; notice that there is 
no contradiction with the no-hair theorem because the Einstein equations are modified by 
the presence of (5i/2- 

In the case 7 > — 1 the total energy is infinite and this may indicate only that solutions 
will not be spherically symmetric. For example the solution may be unstable under axially 
symmetric perturbations and drop to a flux tube in a sort of mass confinement scenario. 

Indeed, regarding stability, even for the physical case 7 < —1 the final word would be 
given by studying the small fluctuations also around the exact solution, to check that the 
non-propagation of the (ghost) scalars and vectors is preserved on a nontrivial background. 
To this aim, we have discussed how the spontaneous Lorentz-breaking is present also in 
the nontrivial background, where we note that in general also rotations are locally broken. 
One expects this also to be a generic feature of massive gravity. 

We showed that we can reach the LI phase by tuning — > 1 in the exact solutions, 
and this results into a well behaved phase, though not the Fierz-Pauli one (gravitons are 
massless) . The fate of the discontinuity and Vainshtein claim for the PF case is thus still an 
open problem and exact solutions of type II with c = 1 are presently under investigation. 
Finally, it would be interesting to speculate on the role of the mass screening in cosmology, 
that may change the form of the Hubble expansion. 
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A. Background for solvable potentials 

For the solvable potential (p.^), we report here the biflat solution as it results from solving 
the EOM ( p. 12 ), as well as the and A2 constants of the linearized analysis. 

The fine tuning conditions to ensure flatness /Ci = /C2 = turn in two relations for the 
two cosmological constants. Defining a„ = (^"^"a^, bn = ui'^^bn and a„ = {dn — (?hn){(? — 
l)/c^, and (3n = {dn + c^bn), we have: 



3/Ci = - 2cj-2^i 



-3/2 



(C2 - 1)7 



+ C 



-3/2 



3/C2 



-1 



8c2(3c2 + 7 + l)a2 - 12c^((7 + 6)0^ + 37 + 2)03 
(6c^ -2)P2 + 24(c2 - - doc^ + 364 - 5a4 

C-5/2 

8((7 + ly + 3)a2C^ + 12((37 + 2)c^ + 7 + 6)030' 



(c^ — 1)7^ 



-5/2 



2{c^ - 3)/?2 + 24(c2 - + doc^ + bbic^ - Sd^) 



(A.l) 



and we remind that one of these is a genuine fine tuning to achieve flatness, as is usual in 
General Relativity, while the other is a complicated equation that may be used to find lo. 
We prefer thinking in reverse and consider uj a free parameter determining the right A2. 
Finally, the lorentz breaking speed of light turns out to be 



di + 4a2 + 603 
61 + 462 + 663 



(A.2) 



The relevant quantities entering in the linearized analysis are the graviton mass A2 and 
the n'^ parameter, that have the following expressions: 



(7 - 2)- 
32 



ao 



C2(c2 -1)(7- 2)7 
.2^1 , 2^ 



X2 = G^^^-^(a2 + 3a3) 



7 



5(c2-l)(7- 2)7(04 + 0^64) 
+ 6c^(l + c^)(7 - 2)7/32 + 32c2(l + c^){j - 2)7/^3 
- 4c2(c2 - 1) [(7 + 2)2^2 + (12 + 47 + 772)03]] I (A.3) 

(A.4) 



B. Interior solution 



For generality we report here the case of a star composed of two spherical regions filled 
with incompressible fiuids of kind 1 and 2 extending from the origin to different radii Ri 
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and i?2; of constant densities pi 2 = /^i^^'^i and negligible pressures. In general we have 
two scenarios, for the three different regions: 

a) for i?2 < ^1 we have < r < R2., R2 < r < Ri ov r > Ri; 

b) for i?i < R2 we have < r < Ri, Ri < r < R2 or r > R2- 

We give only the analitic results for J[r], that is the gravitational potential in gi: 

• Starting from the exterior sohitions r > i?i,2, we find a common value for the exterior 
potential J in both scenarios a) and b): 



J{r) = 1 - ^ 
r 



Mi + 



16/^2 {MiRl 



M2Rf) ^ciJ- 



5(7-2)(7 + l) 
(7-4)(7-2)(7+1)(27-1)ac 



+ 



(B.l) 



• The intermediate solutions are different in the two scenarios: 
for a) i.e. R2 < r < Ri 



^ 2,GMi 2,2G^fi'^uj*M2Rl ^1 \} ~ (2-7)(7+i) ) 48Gr^^/^fi^uj^Mi 

R^^ 5r(2-7)(7 + l)K ^ 5(7-4)(7 + 3)i?j 

96G^/c^l^u^MlR'lr^-'y 9QG^n^uj'^M2R\~'^r'^ 



(7 - 2)(7 + 1)(7 + 3)(27 - 1) ^ (7 - 4)(7 - 2)(7 + 1)(27 - l)n 
for 6) i.e. Ri < r < R2 



(B.2) 



r 5r(2-7)(7 + l) (7 - 2)(7 + l)Kii:2 6(7 - 4)(7 + 3)Ki?3 



(B.3) 



(7 - 2)(7 + 1)(7 + 3)(27 - 1)k (7 - 4)(7 - 2)(7 + 1)(27 - 1) " 
• The inner solutions r < Ri^2 have again a common form in both scenarios a) and b): 



R 



ASGy/cti 



Ml 



+ 



1 



(7-2)(7 + l)i?i Rl) (7-2)(7 + l)Ki22 



16^/i2u;4M2 



+ 



5(7-4)(7 + 3)A^i??i?3 ' (^_2)(7 + 1)(7 + 3)(27-1)k 

As one checks, the solution is regular at the origin. 



.(B.4) 



C. Energy integrals 

In the presence of a time-like Killing vector in GR on can define the notion of total gravi- 
tational energy as a flux from an asymptotic 2-surface that involve only the gravitational 
field at large distance, far from the sources. Consider the following metric 

ds^ = -C{r) dt^ + 2D{r) drdt + A{r) dr^ + B{r) dO^ , (C.l) 
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with the time-Uke Kilhng vector K = Prom the Kihing equation we have 

V% = 0, J^ = nK^. (C.2) 
The Komar energy £ is defined by 



£ = w 

it=ti 



f ^n^J^, (C.3) 
Jt=U 



where h^^ is the induced metric in the hyper-surface t = const, with unit normal and 

t/; is a normalization constant. According to Stokes theorem, given a 3-surface V, for any 
antisymmetric tensor F^i, we have 

[ (fxVhni'V''F^^= [ (fx^(n'^v'^-n^v°')F^p, (C.4) 
Jv JdV ^ ' 

where is the unit normal to dV and 7^/3 is the induced metric in dV . Then Stokes 
theorem gives 



£ = ~ [ V7 (n'^v^ - v'^nP) ^aKp . (C.5) 

^ Jt=ti,r=ri. ^ ' 



In general the Komar energy will depend on 2-surface that bounds the i = const, slice. 
Indeed, from Einstein equations it easy to show that the difference A£ between the Komar 
energy computed with two different bounding 2-surfaccs Si and S2 is proportional to the 
integral of the Ricci tensor over the 3- volume bounded by Ei and E2. As a result the 
Komar energy does not depend on S in a region where the Ricci tensor is vanishing. This 
is indeed the case in a region far from any source. Then 

£ = f Vl(n"vl'-v''nf)Vo.Kfs; (C.6) 

Jt=const.,r—nx>. 

for the induced metric on the 3-surface t = const, and its normal n we get 
dp = A{r) dr^ + B{r) dn^ ; 

(C.7) 



n 



= (AC + f-A'/^^ + DA-'/^4-) 

\ at or J 



and for the induced metric on t, r = const, and its normal ( v is normalized with h ) 

ds'^ = Bir) dn\ v = A^/^-^ . (C.8) 

or 

We have then 

f = - lim . (C.9) 

'•-00 ^/wTac 

One can recover the same result using the language of differential forms. Introducing 
the 1-form J = —{Sd + dS)K in terms of the 1-form K associated with the Killing vector 
K. Prom the Killing equation 6K = 0, then 

J = + dS)K = 6dK = *d*dK = UK^dx'' . (C.IO) 
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Now 



and finally 



d*J = 



*J is time-independent 

t=const. 



(C.ll) 



£ = —w I *J 
't=ti 



w I *6dK = w * * d* dK = w d* dK 
t=ti Jt=ti Jt=ti 

w I *dK = j y/g^g''^ g""" e^^^ap dx'^ A dxf^ 

t=ti,r=ri Jt=ti,r=ri ^ 

ilTW BC 



D. Simplest bigravity 



(C.12) 



Here we analyze the system in the simpler particular case when V \s a. function of q only: 
V = f{q). The Bianchi identities (2.8) for Qi^2-, can be written as^^ 



d^V -[d,\ogq + 2{Vi^-V2u)]{V'xY^ = Q (D.l) 
8(Vi, + V2.)(y'X);;-ya^logg = (D.2) 
and because in the case at hand (V'X)'^ = f'qS'^^, the only non-trivial equation is 



16 



d^ 



diloi 



;f-f 



df,q = 0. 



(D.3) 



Thus, either q =const. or y = Vq = Ciq^^^ + C2q~^^^- However, {gig2)^^^yo = Cl^/gT + 
C2\f^ would imply that the two sectors do not see each other and we are left with two 
independent copies of GR -|- cosmological term. The theory with q =const is the simplest 
of all possible bigravity theories, the EOM reduce to 



with 



M2 2i?2^;+/c2<5^; = -^2^ 



/C2 = g-'/' [/(<?) + 4g/'(<z)]. 



(D.4) 
(D.5) 



(D.6) 
(D.7) 



The effective cosmological constants are thus related; moreover, this simplest bigravity, due 
to the constraint q =const, is equivalent to a single GR -|- unimodular GR, and the two 
sectors share the conformal mode. Finally, besides the diagonal diff also two independent 
volume-preserving diffs are present. 



^^And recall that for any vector field v^^ one has (V2n — Vip)f^ = C^^^o-, with the tensor C^^ = 
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As an example of exact solution, we present the solution for the potential 



V = trlnX = IndetX 



(and we may add also the two cosmological constant terms Aig"-*^/^ and h.2q^^^ to achieve 
flatness). With this potential we have V' X = 1 by construction. The solution in general 
is Schwarzschild-deSitter for both metrics, but (72 is in a different gauge: 



J 
C 

B 
A 



Ai 1 



r 



A2 ( 1 - 2— + C2P^ 
P 



K 

P- 



Ai/J, 

(r3 + A3)1/3 



2 2 

UJ p 



free , 



D^ + AC = A2 
4Ao 



{B 



i\2 



c^iu^Ai {p' 



B 



'\2 



UJ 



2Ai 



(D.9) 
(D.IO) 

(D.ll) 
(D.12) 



This is a family of solutions because A{r) is a free function (!), remnant of the spatial diffs. 
The determinant AC + is fixed by B{r), and for A 7^ it is not constant, at finite 
distance. Then one can also use A to set -D = and get a bidiagonal solution like (^). 
Notice that uj"^ and A2/A1 are free constants, and so also the relative speed of light (? is 
free. 
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